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1 INTRODUCTION 



ABSTRACT 

The problem of detecting dark matter filaments in the cosmic web is considered. Weak 
lensing is an ideal probe of dark matter, and therefore forms the basis of particularly 
promising detection methods. We consider and develop a number of weak lensing 
techniques that could be used to detect filaments in individual or stacked cluster 
fields, and apply them to synthetic lensing data sets in the fields of clusters from 
the Millennium Simulation. These techniques are multipole moments of the shear and 
convergence, mass reconstruction, and parameterized fits to filament mass profiles 
using a Markov Chain Monte Carlo approach. In particular, two new filament detection 
techniques are explored (multipole shear filters and Markov Chain Monte Carlo mass 
profile fits), and we outline the quality of data required to be able to identify and 
quantify filament profiles. We also consider the effects of large scale structure on 
filament detection. We conclude that using these techniques, there will be realistic 
prospects of detecting filaments in data from future space-based missions. The methods 
presented in this paper will be of great use in the identification of dark matter filaments 
in future surveys. 

Key words: cosmology:observations - cosmology: theory - galaxies: clusters: general - 
galaxies: photometry - gravitational lensing 



In a universe dominated by Cold Dark Matter (CDM) such 
as our concordance cosmology, ACDM, hierarchical struc- 
ture formation takes place with the growth of collapsed ob- 
jects progressing via merging of smaller objects. Simulations 
of structure formation show that galaxy clusters are located 
at the intersections of filaments, which form a web- like struc- 
ture throughout the universe. The processes of merging and 
continuous accretion of surrounding matter are thought to 
occur highly anisotropically, with the filaments channeling 
matter along preferre d directions. Moder n galaxy surveys 
such as 2dFGRS (e g. [Colless et all hOQlh) and SPS S (e.g. 
lYork et all (|2000l '). lAdelman-McCarthv et al] (|2008h ) have 
provided a dramatic picture of this so-called 'cosmic web', 
where voids lie between filame ntary arms traced by galaxies. 
From a theoretical perspective. lBond" et al.l (|l996i ) developed 
a theory for the cosmic web in a CDM cosmogony, where fil- 
aments are a product of a primordial tidal field evolving 
under non- linear efi'ects. 
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Detection of inter-cluster filaments has important impli- 
cations for cosmology, confirming the prediction that struc- 
ture grows highly anisotropically from small initial pertur- 
bations in the cosmic density field. Their detection provides 
an important validation of our picture of hierarchical struc - 
ture formation and cluster evolution l|Colberg et al.ll2005l ). 
Dark matter filaments also have a significant role to play 
when considering the total mass budget of the universe - 
the galaxies and gas channeled along these filaments may 
account for as much as half of the baryonic mass budget in 
the universe. 

Inter-cluster filaments have been the target of a num- 
ber of observational searches. The earliest of these used 
bremsstrahlung X-ray emission to det ect hot gas bein g chan - 
neled along filaments. For example, iBriel fc Henrvl (|l995l ) 
used X-ray data from the ROSAT All-Sky Survey, but 
failed to find evidence for filaments. Subsequent X-ra y 
search es have also proved inconclusive (|Scharf et al] (|2000l ). 
iKull fc Bohringer, (,1999ihlDurret et al.l (120031 )'). the main is- 
sue being that it is difficult to ascertain whether strong X- 
ray emission is due to matter in a filament or due to the clus- 
ters interacting. Another proposed observational method 
looks for overdensities of galaxies relative to the background. 
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IPimbblet fc Drinkwate"3 HqoI) reported finding a short fil- 
ament b etween galaxy clust ers A1079/1084 using such a 
method. lEbeling et alJ (j2004l ) report the detection of a fila- 
ment leading into MACS J0717.5-I-3745. However, this ap- 
proach is susceptible to large errors and requires spectro- 
scopic data to establish the redshifts of the target galaxies. 
An alternative detection technique uses a 'skeleton', which is 
defined as the subset of critical lines joining the saddle points 
of a field to its maxima while following the gradient's di- 
rection. The skeleton formalism has been used in numerical 
simulations to extract and analyze the filamentary s t ructu re 
of the universe ()Novikov et al.1 |2005l l. lSousbie et all (j2007h ). 

W ork by iHahn et all (|2007l ') and lAragon-Calvo et al.1 
l|2007l ) demonstrates that it is possible to use the Hessian 
of either the potential or density field to provide a dynam- 
ical classifica tion of filamentary st r uctur e. This work was 
extended by iForero-Romero et al.l (|2009| ). The number of 
eigenvalues with magnitude above a certain threshold at 
each grid point determi nes whether that point belongs to a 
void, sheet or filament. iHahn et al. I ([2007 ) investigated the 
effect filaments have on the properties of dark matter halos, 
and found that filaments influence the magnitude and direc- 
tion of halo spin as well as halo shape. There are a number 
of other filament detection methods, not elaborated here, 
which are discussed in Pimbblet ( 200f^ ). 

One of the most promising methods of filament detec- 
tion is weak gravitational lensing. Since gravitational lensing 
is insensitive to whether the matter is luminous or dark, and 
to its dynamical state, it is an ideal probe of dark matter. A 
number of attem pts have been mad e to image filaments us- 
ing weak lensing. ICrav et "all (|2002l ) claimed to detect a fila- 
me nt in the A901/902 sup ercluster, but this was later shown 
by iHevmans et al.l (|2008|) t o be an artifact of the mass re- 



construction. IKaiser et al.l l|l998h reported a detection of a 
filament between two of the galaxy clusters in the superclus- 
ter MS_0302+17, but th ese results could not be reproduced 
bv lGavazzi et al.1 b004). and hence this detection is called 
into question, loietrich et al.l (|2005| ) used a variety of meth- 
ods to detect a filament between A222/223 - however, they 
point out that it is difficult to objectively define whether 
this lensing signal is caused by a real filament. The Hub- 
ble Space Telescope Cosmic Evolution survey (COSMOS) 
has provide d corroboratory evid ence for a cosmic network 
of filaments (|Massev et al.li2007l '). 

Although there have been a number of observational 
attempts to detect the filamentary cosmic web with various 
degrees of success, conclusive evidence has proved elusive 
in all but a few cases, and widely applicable methods are 
still lacking. In this paper we assess existing techniques to 
detect filaments in weak lensing data and present two new, 
more effective, detection methods. In particular we present 
a method based on multipole filtering of the shear field, and 
a method based on Markov Chain Monte Carlo techniques 
that allows both detection and quantification of filament 
profiles. Motivated by the promise of weak lensing data from 
many thousands of galaxy clusters from large surveys such as 
the Dark Energy Survey (DES) and Large Synoptic Survey 
Telescope (LSST) we also consider detection in stacked data 
sets from a number of clusters. These methods are then ap- 
plied to single and stacked synthetic weak lensing data sets 
from the fields of clusters from the Millennium Simulation 



(jSpringel et al.ll2005h . We will show that weak lensing is an 
ideal and versatile tool to detect dark matter filaments. 

This paper is organized as follows. After a review of the 
basic lensing formalism (section 2), in section 3 we discuss 
the various methods by which we may detect filaments using 
weak lensing. These methods are then applied to data sets 
from the fields of simulated clusters. In section 4 we describe 
these simulations and in section 5 we detail our results. In 
section 6 we analyze the effects of large scale structure noise 
on our results. Finally, in section 7 we present our discussions 
and conclusions. 



2 SUMMARY OF WEAK LENSING 
FORMALISM 

In this section we present the basic formalism of weak grav- 
itational lensing. The surface mass density of a lens at posi- 
tion is denoted by E (G). The critical surface mass density 
of a lens at redshift Zd, with background sources at redshift 
z, is defined as Ecrit = 4Ig"d7d7' where D^, Dd, Dds are 
the observer-source, observer-lens and lens-source angular 
diameter distances respectively. The dimensionless surface 
mass density of the lens is called the convergence, and is de- 
fined by K.{0,z) = . The convergence is related to the 



defiection potential (9) through a Poisson equation, 

2k . (1) 



Whereas the convergence magnifies background objects such 
as galaxies, the shear stretches them tangentially. The com- 
plex shear depends on the second derivatives of the lensing 
potential. 



7 = 71 -I- 272 = (V",!! - V'.aa) /2 -I- it/',: 



(2) 



where commas denote partial differentiation. 

The complex reduced shear is obtained from the shear 
and convergence through the equation g = — n). The 
magnification is given by the inverse Jacobian determinant 
of the lens equation, jJ.{0) — [det A {0)]~^ . Evaluating this 
we obtain. 



(3) 



(1 



We can define the ellipticity of a galaxy, describing its shape 
and orientation, by a complex number with a modulus re- 
lated to the axis ratio r of (1 — r)/(l + r), and phase 
being twice the position angle. In the non-critical regime 
(det ^ > 0) the weakly lensed image of a distant galaxy has 
a complex ellipticity given by. 



(4) 



where e and are the lensed and unlensed complex ellip- 
ticities respectively, and g* is the complex conjugate of the 
reduced shear. Galaxy ellipticities are taken to have a Gaus- 
sian probability density function. 



exp I 



3|2 / 2 



(5) 



Throughout this work we take CTe — 0.2. In weak lensing, 
since the distortion of any individual background galaxy 
is too small to detect, we must look at a large sample of 
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galaxies and statistically detect the presence of dark mat- 
ter, which will manifest itself in, for example, a non-zero 
mean ellipticity of background galaxies. It has been shown 
that the expectatio n of the lensed ellipticity is (e) — g in the 
non-critical regime (ISchramm fc Kavsei 1995) and (e) = 
in the critical regime ifseitz fc Schneider Il995ll . 



3 METHODS TO DETECT FILAMENTS 

In this section we outline the methods of filament detection 
used in this paper: non-parametric mass-reconstruction, pa- 
rameterized model fits using MCMC techniques, and aper- 
ture multipole moments. 



3.3 Aperture Multipole Moments 

ISchneider fc BartelmannI |l993) introduced the concept of 
an aperture multipole moment, defined by. 



d'6» e"u{\e\)e"'"^ti{0o + e) 



(8) 



where U {\0\) is a radially symmetric weight function, and 
we are using polar coordinates {9, <j))- In particular, the aper- 
tu re quadrupole momen t of the convergence has been used 
bv lDietrich et all l|2005l ') as a method for detecting filaments. 
We can express the aperture moments i n terms of shear es- 
timates (jSchneider fc Bartelmannlll997l ), 



3.1 Non- parametric Mass Reconstruction 

We perform finite- field non-p arametric mass reco nstruc- 
tions using the algorithm of .Seitz fc Schneider! l|l996t ). Non- 
parametric mass reconstructions use a statistical analysis of 
many weakly lensed galaxies to identify over-densities in the 
mass distribution. 

Many previous attempts at detecting filaments using 
gravitational lensing have focused on cluster pairs, and in 
this paper we too will consider such filaments. This is be- 
cause filaments between closely separated clusters tend to 
be both strong and regular. However, mass reconstructions 
must be applied with care to closely separated pairs of clus- 
ters. The smoothing length of the mass reconstruction repre- 
sents the typical scale over which the resulting reconstructed 
mass is 'smeared'. If this smoothing scale is greater than 
the inter-cluster separation, this smearing will join the two 
clusters, resulting in a possible false identification of a fila- 
ment. Indeed, this was suspected to have been a factor in 
the false detection of a filament in the A901/ 902 superclus- 
ter jGrav et all l|2002l ): iHevmans et al.l l|2008l )'). 

3.2 Parameterized Mass Models 

The probability distribution for the observed galaxy elliptic- 
ities, pe, can be obtained from the tntrinisc ellipticity dis- 
tribution Pes ([5| as follows. 



pMo) =P^'' (e' (els)) 



g2 



= P.'ie{e\g))- 



* 1 |4 

eg* - 1 



(6) 



Using this, we can compute the log-likelihood function from 
the probability density for each lensed galaxy, pe (e^), 

£-, ^ -J2^np4e4g{e,)) , (7) 

i = l 

which can be evaluated numerically for a trial mass model 
(with corresponding reduced shear field g{0)) given the 
lensed ellipticities. 

For a mass model with a small number of parameters to 
obtain the best-fit parameters it is often appropriate to use a 
minimiza tion technique such as the conjugate gradient algo- 
rtithm of lPowelll (|l964l ). In this paper we instead employ an 
MCMC approach, impleme nting our likelihood c alculation 
as a module in COSMOMC (|Lewis fc Bridle! |2002| ). using it 
as a generic sampler to explore paramete r space with the 
stand a rd Metropolis-H astings algorithm (|Metropolis et af] 
(|l953! ): lHasti"ng3 (|l970l )). 



1 ^ 
(Bo) = iVe"'-* 



U(0i)£ti+1 £xi f , 



(9) 



where n is the number density of galaxies in the aperture, 
and eti = — 7?.(eie~^"*') and e^i = — X(£ie~^"^') are the tan- 
gential and cross components of the shear respectively. U'{9) 
is the derivative of the weight function with respect to the 
radial coordinate. [Dietrich et al.! l|2005!) used a weight func- 
tion of the form. 



u{e) 



1- 







for ^ 9max 



otherwise 



(10) 



where Smax is the aperture radius. This serves well as a gen- 
eral multi-purpose filter, but may not be optimal for filament 
detection, as it does not match their expected mass profiles 
dColberg et al.ll2005! ). 

In reality, a single order of multipole will not completely 
characterize a mass distribution - any given mass distribu- 
tion will in general give non-zero results for a variety of 
different multipole orders. Therefore, instead of applying a 
quadrupole filter alone to detect a filament, a more optimal 
method would use a superposition of filters thus. 



(11) 



where the constants at can be chosen to maximize the signal- 
to-noise. 

We can also define multipole moments of the shear field, 



C'"'(0o)=/ d'0ri7(|6l|)e'"'^7(0o + 6>) . (12) 
Jo 

This can be related to the convergence through the re- 
lation, 



where, 



7 (6») = i j d^e e'v{e - e V(f ) 



^^^> = — w — 



(13) 



(14) 



The integral eq. (|12|) can also be written in terms of 
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shear estimates to give, 
1 ^ 

C'"'(0o) = =ye<"+'''^'ert/(eOh+iex] . (15) 
n ^ — ' 

z — 1 

To calculate the noise associated with either the shear 
or convergence multipo l e filte rs, we use the result of 
ISchneider fc Bartelmannl (|l997l ). Both the shear and mul- 
tipole filters have the form, 



{Go) = / d^e e' 



ni<p 



— ^7t(6»; 6»o) + 1— ^7x (O; 0o] 



(16) 



where bt{9) and bx{0) are functions of the radial distance 
coordinate and 7t and 7x are defined by the following, 



7t(6»;6»o) = -7^ [7 (6» + 6»o) e-"'^ 



(17) 



and 



7x (0;0o) = -X [7(0 + 00) e""''^] . (18) 

Using the result of lSchneider &: Bartelmannl l|l997l ') for 
filters of this form, the corresponding ensemble-averaged 
noise is, 



2 r9 
n Jo 



dy 



bUy) + bl{y) 



y 



(19) 



In this paper, we define the signal-to-noise ratio of a 
multipole filter by. 



(20) 



3.4 Stacking 

Since filaments are very weak, their detection will usually 
require high background galaxy counts, above that typically 
achievable at present. One solution to this problem would 
be to stack catalogues of background galaxies in the clus- 
ter fields. Once these lensed catalogues have been stacked, 
one can apply the techniques described above to search for 
structures of interest. 

In this paper, we use stacked catalogues to increase the 
ease with which we are able to detect filaments. For stacking 
to improve our chances of detection, we need to align fila- 
ments in the stacked data. Since we use stacked catalogues 
to detect filaments that cannot be detected on single cluster 
images, we must infer the probable filament direction be- 
fore stacking the data. There are two possible ways of doing 
this. The first method uses the result that a filament is of- 
tcn aligned with the major axis of a cluster ( IColberg et al.l 
[2005. ) . Since we can identify the cluster major axis from a 
meiss reconstruction, we can then align these major axes 
in our stacked catalogues. Depending on the quality of the 
mass reconstruction or parameterized model fit, one could 
also align cluster centres using optical images (assuming that 
light traces mass). The second method involves using clus- 
ter pairs as 'markers' for the starts and ends of filaments. 
It is expected that cluster pairs will be connected by fila- 
ments, so by aligning the axis running through the cluster 
centres we also align any filament between them. The cluster 



pairs must be close enough to ensure that a filament con- 
nects them, but separated by several virial radii to minimize 
contamination by a cluster signal. 

Both methods have advantages and disadvantages, and 
the most useful method will depend on the data-set being 
analyzed. The first method can be applied to all clusters in a 
sample, whereas the second requires cluster pairs with a par- 
ticular spacing, which are rarer. However, the second method 
has the advantage that filaments between cluster pairs tend 
to be straighter and more regular th an those filaments wh ich 
do not connect to another cluster (|Colberg et al.|[2005l ). If 
we simply stack the major axes, a significant fraction of the 
filaments will be very diffuse, 'warped' or irregular - stack- 
ing such filaments does not increase the ease of detection, 
as there will be little overlap of irregular filaments from 
each cluster. Clusters with similar separations could also 
be stacked to maximize the extent of filament available for 
analysis. In this way, the dependence of filament properties 
on inter-cluster separation can also be assessed. 



4 NUMERICAL SIMULATIONS: CLUSTERS 
AND SOURCES 

For this study we m ake use of the Millennium Simulation 
(|Springel et al.ll2005h . a large cosmological N-body simula- 
tion that follows 2160^ dark matter particles from z = 127 
to 2; = in a periodic box of 500 Mpc on a side. The 
cosmology adopted for the Millennium Simulation (MS) is a 
flat ACDM model with h — 0.73, = 0.25, and a power 
spectrum normalisation on a scale of 8/i~^Mpc (i.e., the rms 
linear mass fluctuation in a sphere of radius 8h~^ Mpc ex- 
trapolated to present-day) of as = 0.9. These parameters 
are consistent with the latest measurements of tempera- 
ture and polarization anisotropies in the cosmic microwave 
background (Komatsu et al. 2009) , although the value of erg 
adopted for the MS is larger than the maximum-likelihood 
CMB value by 2a. The larger value of as means that the 
MS will have more massive clusters than than a universe 
with the WMAP 5-year cosmology. However, for the pur- 
poses of our study this is not important as we are interested 
in the lensing signal (and structure) of individual clusters 
and their fllaments and not the statistics of cluster popula- 
tion as a whole. The mass structure of individual clusters 
(e.g. halo concentration), and presumably their filaments, 
depends only weakly on cosmological parameters (e.g. Duffy 
et al. 2009) and therefore the exact choice of parameters is 
not expected to inffuence our results or conclusions. 

From the MS, we select the 100 most massive galaxy 
clusters at Zc = 0.2 for analysis. These are the systems 
which are expected to have the most prominent filaments. 
The clusters are identified using a standard friends-of-friends 
(FOF) algorithm with a linking length of 0.2 times the 
mean interparticle separation. The clusters we have se- 
lected range in mass from M200 = 3.90 x 10^* Mq to 
M200 ~ 2.16 X 10^^ Mq (for comparison, the mass reso- 
lution of the MS is 8.6 x 10*/i"^Mq.) We extract aU of the 
dark matter particles in a 50/i~^ Mpc box centred on the 
FOF halo centre of mass. For each box, we generate surface 
mass density maps by projecting along the axes of the box; 
we make 3 separate projected maps, by projecting along the 
X-, y-, and z-axes of the box. The dark matter particles are 
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Figure 1. In this figure, the left-hand panels show the simulated cluster convergence maps on a logarithmic scale, whereas the right- 
hand panels show the corresponding mass reconstructions. The background galaxy density is 100 arcmin"^, and the smoothing scale 
is 1.5/i^^Mpc. The side length of the field is 3.5°, and the reconstruction is calculated on a 100 X 100 pixel grid. We see that the 
reconstruction provides some evidence for strong filaments, but is unable to detect weaker filaments. 



interpolated onto a regular 2D grid (map) using a triangular 
shaped clouds interpolation alogrithm. We fixed the num- 
ber of pixels in the 2D grids to be 4000^^, implying a pixel 
length of 12.5 kpc, which is comparable to the gravita- 
tional softening length of the MS (which is 5/i~^kpc). As 
we shall see, this resolution is more than sufficient for mak- 
ing weak lensing predictions of clusters and their filaments 
for current and future (e.g., space-based) observations; the 
number densities of sources imply effective resolution well in 
excess of several tens of kpc. 

Simulated galaxy catalogues are created by assigning 
galaxies random positions in the field, and by simulating 
galajcy ellipticities as Gaussian random deviates. A magnifi- 
cation bia s is als o included in the simulated population (e.g. 
ICanizarei (|l982l )l. which changes the galaxy number counts 
after weak lensing. The local cumulative number counts 
n {G; S) above a fiux limit S, are related to the unlensed 
counts no (S) by. 



then. 



n{e-S) = -no (- 



(21) 



where n is the magnification. If we assume that the number 
counts locally follow a power law of the form no oc S~^, 



n (0) = nop (0) 



0-1 



(22) 



at any fixed threshold. This implies that if the intrinsic 
counts are flatter than 1, then the lensed counts will be re- 
duced relative to the unlensed counts. In this paper we use 
P = 0.5, typical of the faint galaxy populations studied in 
weak lensing. To account for the effect of magnification we 
reject a given galaxy from the catalogue if a uniform ran- 
dom deviate in the range [0, 1] is greater than ■ This 
means that the higher the magnification of a given galaxy, 
the more likely it is to be excluded from the catalogue. 

In our analysis, we assume that the filaments associ- 
ated with a cluster are at the same redshift as the cluster 
itself. Thus our typical simulation box of 50 Mpc is pro- 
jected onto a lens plane at the cluster redshift, — 0.2. The 
cluster resdhift is sufficiently low such that the redshift dis- 
tribu tion of the source galaxy population can be neglected 
(e.g. iKing fc Schneider] (|200ll )1: we take the source galaxy 
redshift to be 2g = 1, characteristic of fairly deep ground- 
based observations. 
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5 RESULTS 

We investigated several metliods for filament detection: mass 
reconstructions, multipole filters in both the convergence 
and shear, as well as MCMC fits to filament profiles. In 
this section we describe the results of these techniques. In 
each case, we analyzed a square field centred on the target 
cluster, with a side-length of 3.5°. 



5.1 Mass Reconstruction and Aperture 
Convergence Multipole Moments 

We performed finite field mass reconstructions on single 
galaxy cluster synthetic data sets using background galaxy 
densities of 30 arcmin"^ and 100 arcmin"^, the former being 
typical of current ground-based observations and the latter 
being representative of future space-based missions. The re- 
sults are shown in Fig. 1. With a background galaxy density 
of 30 arcmin"^, filaments are not detected in reconstruc- 
tions. With a background galaxy density of 100 arcmin"^, 
strong filaments are detected more clearly, but weaker fil- 
aments remain undetected. These strong filaments are all 
characterized by having regions of the filament which con- 
tain large clumps of dark matter (corresponding to massive 
galaxy groups or poor clusters), and it is these regions that 
the mass reconstruction is able to detect. In Fig. [T] there is 
visual evidence for particularly strong filaments in clusters, 
but not for weaker ones which are missed. 

The multipole filter outlined in eq. (9) was applied to 
the same data sets, using n = 1 and n = 2 (a dipole mo- 
ment and a quadrupole moment). A variety of filter func- 
tions and aperture sizes were tried, but no filaments were 
detected by either multipole moment. We also superposed 
the two multipole moments as per eq. (|ll|l , but were unable 
to make a filament detection. We lowered the ellipticity dis- 
persion in the synthetic galaxy catalogues to — 0.05, at 
which point both Q^^^ and Q^^^ registered weak detections 
of the strong filaments. Although this CTe is unrealisitic, it 
allows us to make a rough estimate of how many clusters 
we would need to stack for these multipole moments to be- 
come effective. Taking the noise to vary as where A'^ is 
the number of stacked clusters, we estimate that we would 
need to stack approximately 20 clusters to detect filaments 
using this method. This is likely to be a lower bound, due to 
imperfect filament alignment and the presence of irregular 
filaments in the stacked clusters. 

The mass reconstruction was then applied to stacked 
data fields using both of the stacking methods outlined 
in section 3.4. As predicted, stacking using double clusters 
proved to be more effective, and in Fig. [2] we show the mass 
reconstruction for 10 stacked galaxy fields, rotated such that 
the axis joining the centres of the two clusters lies in the 
vertical direction. To select cluster pairs, we identified all 
those fields for which one cluster was observed to be within 
25/i~^Mpc of another cluster. The density of background 
galaxies used for each cluster was 30 arcmin"'^, and the 
smoothing scale used was 1.5/i~'^Mpc, chosen to approx- 
imately correspond to the predicted width of a filament 
l|Colberg et al.1 (|2005^ '). We see clear evidence of filamentary 
extension in the inter-cluster region. The region shown in 
Fig. [2] lies outside the virial radii of both clusters - in other 
words, clusters don't contribute to the signal. A mass re- 




Figure 2. Mass reconstruction for 10 stacked clusters with a 
background galaxy density of 30 arcmin"^. Each stacked data 
set corresponds to a cluster pair, with each individual data set ro- 
tated such that the cluster-cluster axes are aligned in the stacked 
catalogue. The smoothing scale used in the reconstruction was 
1.5 Mpc. The region shown above is the inter-cluster region, 
with a vertical height of 1.3° and a horizontal length of 3.5°. We 
see clear evidence for filamentary presence between the cluster 
pairs. 

construction could also detect evidence for filaments when 
major axes of clusters were stacked in the vertical direction, 
but this signal was weaker than that shown in Fig. [5] 

We applied the first and second multipole moments Q'^-* 
and Q'^' to the stacked data. We detected a faint Q'^-* signal 
and no Q'^' signal. 

5.2 Shear Multipole Moments 

The shear filter fea. ll5|l was applied to a sample of clusters 
using the n = multipole and background galaxy densities 
of 30 arcmin"^ and 100 arcmin"^. The signal to noise maps 
corresponding to two cluster fields and an aperture radius 
of 2.5 h^^ Mpc are presented in Fig.O As before, the aper- 
ture mass was chosen to approximately correspond to the 
expected width of a filament. 

The filtered results, in Fig.O show clear evidence of fil- 
ament detection. As expected, as galaxy density increases, 
filament detections are made at a higher signal-to-noise, 
and weaker filaments are more apparent. The strong fila- 
ments can be seen at both galaxy densities - at 30 arcmin"^ 
the signal to noise along these strong filaments is approx- 
imately 3-5, and in certain regions of the filament even 
higher; at 100 arcmin"^ the signal-to-noise along the length 
of these strong filaments is in the region 5-7, a very ro- 
bust detection. Weaker filaments are less apparent at lower 
background galaxy densities, although in Fig. [3] there is ev- 
idence of weaker filaments in the 30 arcmin"^ images. At 
100 arcmin"^ even these weak filaments come into view at 
signal-to-noise typically in the range 3-4 along the length of 
the filament. 

The aperture size should be chosen to correspond 
roughly to the scale one wishes to investigate. Aperture sizes 
that are too small will produce weak signals, whereas if the 
aperture is too large it will pick up signal from mass un- 
related to the structure of interest. In the specific case of 
detecting filaments, which are much weaker than the clus- 
ter, any signal sufficiently close to the cluster should not be 
regarded as necessarily corresponding to filament detection, 
as such filters will enclose some of the cluster as well as any 
putative filament. To be confident of filament detection, the 
signal must at least extend outside r2oo + ^max, where r2oo 
is the virial radius of the cluster and rmax is the aperture 



Probing the cosmic web: inter-cluster filament detection using gravitational lensing 7 





— 




■ 1^ 






> 


■■■1 








1 









I 

1 





Figure 3. Results for the n = shear filter using an aperture radius of 2.5 Mpc. The top panels show the cluster convergence 
map, on a logarithmic scale. The middle panels show a signal-to-noise map when the shear filter is applied to the cluster field with a 
lensed galaxy density of 30 arcmin"^. The bottom panels show the same result, but this time for a galaxy density of 100 arcmin"^. The 
contours correspond to signal-to-noise values of 3, 4 and 5. The fields are of side length 3.5°, and the grid size upon which the filtering 
is performed is 100 X 100 pixels. 



size. In Fig. [3] we see that the asserted detections satisfy this 
criterion. 

As discussed in earlier sections, using a single multipole 
is not the optimal method for detection. By superposing 
many multipoles, one can choose the coefficients in the ex- 
pansion to maximize signal-to-noise. On investigation how- 
ever, adding higher multipoles made only a small improve- 
ment to filament detection and is unnecessary. One may also 



hope to improve the significance of a detection by altering 
the filter function, however, in this case attempting to match 
the filter function to expected filament profiles did not lead 
to a great improvement. The weight function in eq. (10) is 
suitable for filament detection. 
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5.3 Application of MCMC Techniques 

We now consider the effectiveness of MCMC techniques in 
detecting and describing filaments using a simple parame- 
terized bar m odel, with a p r imary focus on those between 
cluster pairs. IColberg et atl (|2005l ) identified filaments be- 
tween neighbouring clusters in the N-body simulations of 



iKauffmann et al.l l|l999l '). They found that approximately 
40% of all filaments between clusters are straight and on 
centre with respect to the cluster-cluster axes; we shall term 
such filaments 'regular'. Other filaments are 'warped', ir- 
regular or off-centre. The former category of filaments are 
amenable to parameterized mass model fits. 

To select a sample of regular filaments, we selected pairs 
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of clusters in our fields separated by a distance of less than 
Clusters such as this can be easily picked up 
on mass reconstructions, even at relatively low lensed back- 
ground galaxy counts. Filaments connecting clusters tend to 
be straighter and more regular compared to those filaments 
which do not connect cluster pairs. A further advantage of 
focussing on filaments connecting 'closely' spaced clusters is 
that since the clusters provide a start and end point to the 
filament, this means we can constrain to a good approxima- 
tion both the filament's orientation and the position of its 
central axis. This reduces the number of parameters we have 

to fit in our mass model. 

As noted above we use COSMOMC (|Lewis fc Bridle! 

l|2002l )') with a Metropolis algorithm sampler. We dynam- 
ically update the proposal density from the progressive co- 
variance of post-burn-in samples, and we use standard con- 
vergence tests. Analysis of the chains produced uses GetDist 
distributed with COSMOMC. Before fitting a mass model to 
the filament itself, NFW profiles were fit to the cluster ha- 
los, yielding values for their virial radii r^ir- We then used 
all galaxies which were a distance greater than two virial 
radii away from the cluster cores. This was to ensure that 
no artificial contribution to our bar fits was made by the 

clus ters themselves . 

IColberg et al.l (|2005l ) investigated the mass distribu- 
tions of the filaments in their simulations, and determined 
that they tend to have cores of approximately constant sur- 
face density, after which their profiles are well approximated 
by a decline. Therefore, we fit the profile, 

«(0 = 7^ , (23) 

where r is the perpendicular distance from the centre of 
the filament, and the amplitude of the model bar, ko, and 
the core radius, rc, are the free parameters. The shear field 
corresponding to this bar can be derived using symmetry 
arguments. If we take one axis as lying along the filament 
(labelled '2') and the other axis orthogonal to the filament's 
length (labelled '1'), we know that the deflection potential 
cannot vary with respect to the former coordinate. This 
means we can set all derivatives with respect to this axis 
to zero in eq. (2). It is then clear that in this case, 

71 = «; 72 = . (24) 

We use flat priors on parameters, covering as broad a range 
as is physically possible, errors also being accounted for. 
In Fig. m we show the likelihood contours obtained for three 
clusters, using a background galaxy density of 100 arcmin"'^. 
A fit to a 'blank' field (a region of the simulation associated 
with neither the cluster nor a filament) was also performed, 
and as expected the MCMC fi t gave core rad ii and bar am- 
plitudes consistent with zero. IColberg et all (|2005l ) report 
that for straight filaments, the typical core radius of a fila- 
ment lies between 1.5 and 2.0 Mpc and our MCMC fits 
are consistent with this result. 



^ Note that in the z ^ 0.2 snapshot of the Millennium simula- 
tion, almost 25% of galaxy clusters with virial mass in excess of 
2 X lO^'*h~^M0 (mass inside a radius where the mean density 
exceeds 200 times the critical density for the simulation) have a 
neighbour exceeding this mass limit. 



The results were repeated for lower galaxy densities (30 
arcmin"^), and as expected, the general effect of this is to 
increase the spread in our parameter distributions. Although 
the fits were still inconsistent with an empty field, in the 
majority of cases the parameter range was so broad as to 
not provide physically meaningful constraints. 

The filaments analyzed in Fig. U all connect cluster 
pairs, but we may also apply this MCMC technique to fila- 
ments that do not connect closely spaced clusters. We took 
clusters which were not part of a cluster pair and applied 
the shear filter and mass reconstruction to the data fields, 
using a background galaxy population of 30 arcmin"'^. We 
then identified filament candidates from the results, and fit 
filament profiles to these candidates. We found that the re- 
sulting profiles gave core radii and amplitude in the expected 
range, with neither quantity consistent with zero. 



6 THE EFFECTS OF LARGE SCALE 
STRUCTURE 

As weak lensing is sensitive to all matter along the line of 
sight, the large scale structure (LSS) will make a contribu - 
tion to any lensing signal. For instance, iHoekstral (|2002| ) 
demonstrated that LSS is a major source of uncertainty 
in determining the virial mass and concentration param- 
eters of clusters. This was confirmed by iDodelsonI (|2004 ). 
who also presented techniques to minimize the effects of the 
LSS. Given this previous work, we expect the presence of 
LSS to have a deleterious effect on the results presented in 
this paper. To investigate the influence of LSS on our re- 
sults, we simulated Gaussian LSS fields with a convergence 
power spectrum corresponding to the cosmology used in pre- 
vious sections, and added them to our cluster simulations. 
The convergence rang e of these LSS fi elds are consistent 
with the findings of White fc Vald (|2004l ). excluding regions 
where there are non-linear structures. The inclusion of non- 
linear structures will thus increase the scatter beyond that 
of the analysis below. Additionally, large structures located 
along the line of sight could bias mass estimates. However, 
with photometric redshift estimates for sources in the field 
one could identify spikes in redshift space that are likely to 
impact on the results, and partially account for them in the 
models. 

To quantify the eff ects of LSS on the shear filters, we can 
extend the results of iSchneider et al (|l998:l. In the presence 
of LSS, the ensemble average of the shear monopole will be 
unchanged. The dispersion due to LSS can be calculated as 
follows, 

((C'°' (0))') = / d^0'U{0') I d^0{j{0'h{0)) , (25) 

where the shear correlation function c an be related to 
the p ower spectrum as described in e.g. ISchneider et al.l 
i2003), 

(7(0)7(0))= P ^Mie)[PE{l) - Pb{1) +2iPEB{l)] . 

(26) 

In this equation, Pe and Pb are the E- and B-mode 
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power spectra outlined in ISchneider et al.l l|2002h . The dis- 
persion due to different orders of multipole can be calculated 
in an analogous fashion. 

When the shear filter was applied to the clusters with 
the addition of LSS noise, it was found that this noise con- 
fused and obscured filament detection. In particular, at a 
galaxy density of 30 arcmin"^ filaments are difficult to dis- 
tinguish from spurious signals due to the LSS noise - at 
this galaxy density we were unable to make any clear detec- 
tions. The problems from the LSS are mitigated somewhat 
at higher galaxy densities, with strong filaments (i.e. those 
between cluster pairs) detected clearly at a galaxy density 
of 100 arcmin"^. 

To investigate the effects of LSS on our MCMC anal- 
ysis, we simulated a number of LSS noise realizations and 
used MCMC to find the best-fit parameter (the values for 
which there is a peak in the marginalized ID distribution). It 
was found that the addition of LSS noise caused a spread in 
the values of the best-fit parameters, but there was no bias 
(i.e. the best-fit parameters were not systematically higher 
or lower) . For instance, repeating an MCMC analysis of the 
cluster considered in the top left hand panel of Fig. 4 seven 
times we found that the width ranged between 2.8-5 Mpc, 
and the amplitude ranged between 4.5-6.1 xlO~^. The re- 
sults are still consistent with the presence of a bar; however, 
if the results are to be used to quantitatively analyze fila- 
ment profiles the dispersion in the results needs to be re- 
duced. In practice this could be done by statistically model- 
ing the LSS structure noise and including it in the likelihood 
function. 



7 DISCUSSION AND CONCLUSIONS 

The typical strength of an intercluster filament is k ~ 0.01, 
compared to a convergence in the central cluster regions of 
K ~ 0.5. Therefore, the reason filaments are so difficult to de- 
tect using any method is due to the fact that they are much 
weaker than the cluster itself, and are quickly 'lost' in the 
noise originating from the intrinsic ellipticity dispersion of 
the background galaxies. In this paper, we have shown sev- 
eral methods that can be used to reliably detect filaments. 
By implementing a variety of methods, we have also shown 
how versatile weak gravitational lensing is as a tool for de- 
tecting filaments - each method described in the paper has 
different strengths and weaknesses that make them suited 
to different data sets. 

The mass reconstruction is one of the most common 
methods to detect overdensities in a mass distribution. 
W e applied the finite field mass reconstruction algorithm 
of ISeitz fc Schneiderl (|l996h to both single and stacked 
clusters. For single clusters, reconstructions using a back- 
ground galaxy density of 30 arcmin"^ weakly detect only 
the strongest filaments, with the majority of filaments unde- 
tected. As background galaxy density increases, so does the 
ease of filament detection - at 100 arcmin"^ we are able to 
identify some filaments in reconstructions. The mass recon- 
struction is most effective when used with stacked data sets. 
The stacking method has an impact on the quality of the 
results, with the double-cluster method proving more effec- 
tive than aligning cluster major axes. A major disadvantage 
of mass reconstructions is that they can give false filamen- 



tary structures due to the smearing effect of the smoothing 
scale. This is not a problem unique to mass reconstructions, 
and any method which smoothes results over certain scales 
in an aperture will suffer from the same difficulty. This is a 
particular problem for closely-spaced clusters, which will ap- 
pear to merge in the reconstructed image if the smoothing 
scale is greater than the inter-cluster separation. The sec- 
ond major disadvantage of a mass reconstruction is that it 
does not allow easy calculation of the significance of any fil- 
ament detection. The origin of this effect lies in the fact that 
in any reconstructi on, the shear field must be smoothed to 
avoid infinite noise ijKaiser et al.lll993 l. This leads to corre- 
lated errors in the resulting convergence map, making error 
bars very difficult to construct. Thirdly, because filaments 
are weak, they are often located in reg ions that are domi- 
nated by noise. As lDietrich et all (|2005l l points out, in such 
a scenario the value of the mass sheet degeneracy can fiuc- 
tuate by as much as the filament strength between different 
galaxy populations, making any filament detection in an in- 
dividual cluster more questionable. This third point does 
not apply for the stacked clusters, as the stacking ensures 
that noise is reduced to sub-dominant levels even in fila- 
mentary regions. Therefore, it is advisable to use the mass 
reconstruction in conjunction with other methods described 
in this paper, both to confirm any detection and to allow 
the significance of any detection to be accurately assessed. 

In general, the multipole moments of the convergence 
peformed poorly, and were unable to register any strong fil- 
ament detections in single clusters. Although IPietrich et al.l 
(2005) used Q'-^-' to provide evidence for a filament in 
A222/223, this cluster pair was very closely spaced and the 
corresponding filament was thus much stronger than any 
considered in this paper. It is also unclear how much of a 
contribution the filament makes to the quadrupole signal, 
as it is likely that the majority of the signal comes from the 
clusters themselves. The quadrupole signal fails to detect 
filaments in the stacked data due to misalignments. Mis- 
alignment between filaments will create a 'mass-sheet' in 
the resulting stacked data, which will give zero quadrupole 
signal. 

It is not only possible to define multipole moments of 
the convergence, but also of the shear. Here we have shown 
that the shear monopole is very effective in picking out fila- 
ments. The reason for this is that the shear monopole per- 
forms a (weighted) sum of all shears in an aperture. If these 
shears are coherent over the aperture, then the net result 
will be a large signal. The shear signals corresponding to 
a filament will all point in approximately the same direc- 
tion along the filament's length, and thus filaments are well 
suited to detection using this method. However, this method 
is not tailored to exclusively find filaments, and a variety of 
mass distributions will give a strong signal. The LSS adds 
noise to the shear monopole signal which can, in the cases 
of weak filaments or low galaxy densities, cause confusion, 
false-identification and obscuration of the filament signal. 
However, we have shown that at galaxy densities character- 
istic of future space-based missions, or if particularly strong 
filaments are targeted, the filament signal will still stand out 
clearly against the LSS noise. 

We demonstrated that one can use MCMC techniques 
to fit filament profiles, even at relatively low galaxy den- 
sities. This is an important result as it shows that weak 
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lensing can be used not only for filament detection, but also 
for quantitative analysis of filament properties. In the fu- 
ture, this will allow us to directly compare the predictions 
of simulations to observations. The scatter in the best fit 
parameters caused by LSS has been quantified, and hinders 
the ability of this method to quantify profiles accurately - in 
future work, the LSS should be incorporated into the like- 
lihood function to minimize the i mpact on the result s (e.g. 
as per the prescription laid out in lHu fc White! l|200lh ). The 
MCMC method was also shown to be quite versatile. It can 
be used on its own, as it was when we looked for filaments 
between cluster pairs, or it can be used in conjunction with 
other techniques such as the mass reconstruction or in par- 
ticular the shear filter to first of all search for a filament can- 
didate. There are a number of caveats associated with using 
MCMC fits. Firstly, it is model-dependent - this involves 
assuming a pr ofile for the filament, in this paper based on 
the findings of IColberg et al.l ()2005h . Models that more ac- 
curately refiect the 'true' filament profiles will achieve better 
likelihood values. Secondly, since it is very difficult to write 
down a model for irregular or warped filaments, this method 
is only suitable for investigating the properties of straight, 
regular filaments. This represents only a fraction of the to- 
tal filament population. Thirdly, the method works best if 
we have some prior information on the location of the fil- 
ament, either from one of the other detection methods, or 
by using cluster pairs between which filaments are thought 
to be constrained. Without a prior on the filament position, 
fitting any profile would become a complicated procedure. 

Stacking is a technique that can be used in conjunction 
with any of the methods described above, to improve the 
chance of filament detection. Stacking is particularly use- 
ful when each individual cluster lenses a low background 
galaxy density, but the survey contains many such clusters. 
In this paper we suggested two possible stacking methods: 
the first aligned cluster major axes in the stacked images 
and the second aligned axes joining cluster pairs. The lat- 
ter method is far more effective. This is because, although 
filaments are often observed to align with the major axis of 
a cluster, there is no guarantee that they are either regular 
or strong. In almost all cases there is imperfect alignment 
of the filaments with the cluster major axis, and sometimes 
there is no filament aligned in this direction at all. If irreg- 
ular or misaligned filaments are stacked the net effect will 
be to produce a low density mass sheet. This is not ideal for 
any of the methods described in this paper. On the other 
hand, filaments between cluster pairs tend to be straight, 
strong and regular and hence are easily stacked. The only 
disadvantage of using cluster pairs is their relative scarcity 
- whereas any cluster can be used when aligning the major 
axes, fewer cluster pairs will exist in a survey. 

In summary, in this paper we have presented a vari- 
ety of new detection methods, some of which (for example 
the shear filter and MCMC) give promising results. How- 
ever, due to the low surface densities of filaments, and the 
additional effects of LSS noise, the task of detecting fila- 
ments in actual data remains difficult. In reality, filament 
detection will either require the high galaxy densities of fu- 
ture space-based missions, or the stacking of a number of 
clusters. Current space- based lensing studie s achieve « 80 
galaxies arcmin"^ (e.g. l|Leonard et al.|[2007t )). and we could 
achieve in excess of 250 arcmin"^ for deep space-based ob- 



servations (e.g. (jRhodes et al.|[2004 )). Although there is a 
trade-off between the area and depth of a survey, so that 
those targeting many thousands of square degrees will typi- 
cally not go so deep, the most massive clusters such as those 
considered in this work will be prime targets for deep tar- 
geted observations. Stacking the data from the fields of less 
massive clusters in non-targeted large area surveys will also 
be extremely useful in statistically constraining the proper- 
ties of filaments. 

There is one further method, not discussed in this pa- 
per, that could prove valuable in detecting filaments - grav- 
itational fiexion. In traditional weak lensing studies, such as 
the formalism used in this paper, the lens equation is ap- 
proximated as linear. This means that any distortion in the 
ellipticity of a galaxy will be aligned purely tangentially to 
the lens. Physically, such a linear approximation is equiva- 
lent to assuming no variation in the lens field over the scale 
of the lensed image. However, if we do not make such an as- 
sumption, the lens equation becomes non-linear, and in addi- 
tion to the convergence and shear, gravitational distortion is 
determined by the first and second fiexion. Taken together, 
fiexion introduces 'arciness' to the ima ge and some radial 
alignm ent with respect to the lens (see IColdberg fc BaconI 
pOOSi ') fo r further details ) . It h as been demonstrated, for ex- 
ample bv lLeonard et af] (|2008l ) , that gravitational flexion is 
an effective probe of substructure in galaxy clusters. This 
suggests that fiexion may also be a useful tool for filament 
detection. 

Besides detecting filaments using weak lensing data, as 
we have noted in the Introduction, there are also a vari- 
ety of other complementary detection methods. Sunyayev- 
Zeldovich observations could in future be used in combina- 
tion with more established methods of filament detection 
(e.g. weak lensing. X-ray, galaxy overdensity) to provide 
more optimal algorithms for detection and even stronger 
constraints on profiles. The area of overlap between DES 
and the South Pole Telescope, SPT, will result in several 
tens of thousands of galaxy clusters in common. The ability 
to probe the total matter content of clusters and their en- 
vironments, coupled with deep optical and infrared imaging 
and spectroscopic observations, will enable us to determine 
galaxy bias as a function of environment. It will also im- 
prove our understanding of the formation of the most mas- 
sive bound objects in the universe. Weak lensing of back- 
ground galaxy populations becomes le ss effective for use as 
a probe of higher redshift structures l|Lewis fc Kinel[2006l ) 
and beyond z ^ 1 an alternative source population will be 
essential - such as the 21cm emission from high redshift 
proto-structures that will be imaged with SKA, or the CMB 
as seen by a future mission. The ability to trace the evolu- 
tion of filaments in the cosmic web will have important con- 
sequences for models of structure formation e.g. how much 
material is bound in filaments as a function of epoch. 
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